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Abstract 

We calculate the finite temperature effective potential of at 
the two loop order of the 2PPI expansion. This expansion contains all 
diagrams which remain connected when two lines meeting at the same 
point are cut and therefore sums systematically the bubble graphs. 
At one loop in the 2PPI expansion, the symmetry restoring phase 
transition is first order. At two loops, we find a second order phase 
transition with mean field critical exponents. 
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1 Introduction 



The X(j)'^ theory at finite temperature is of great interest for the field of phase 
transitions in the early universe and heavy ion collisions. When used as a 
simple model for the Higgs particle in the standard model of electroweak 
interactions, it may allow the study of symmetry breaking phase transitions 
in the early universe. For N=4 scalar fields, it is also a model of chiral 
symmetry breaking in Q.C.D. and hence relevant for the theoretical study 
of heavy ion coUisions. Moreover, this theory is an excellent theoretical 
laboratory where analytic non-perturbative methods can be tested. 

Wc recently [1,2] introduced the 2PPI expansion as an analytic pertur- 
bative variational method in quantum field theory. It is perturbativc in the 
sense that one calculates Feynman diagrams which remain connected when 
two internal fines meeting at the same point are cut (2PPI = two particle 
point irreducible). It is variational because the contribution of 2PPR dia- 
grams is absorbed in a variational effective mass. It can be compared to the 
CJT [3] formalism which is an expansion in 2PI diagrams where the varia- 
tional object is the two particle Greens function. The amount of resummation 
in the 2PPI expansion is less than in the 2PI expansion but because the vari- 
ational object is just a mass parameter, it has the advantage of being much 
more tractable at higher order. The 2PPI expansion has also some other 
advantages when compared with other resummation methods. Besides being 
a systematic expansion, its renormalization is straightforward. In [2], it was 
shown that the 2PPI expansion can be renormalized with the usual countert- 
erms in a mass independent scheme. The simplest approach is to calculate 
the 2PPI diagrams using dimensional regularization and subtract the pole 
terms. The 2PPI expansion is especially well suited for finite temperature 
field theory. It does not have the well known problems with perturbation the- 
ory at high temperature because the contribution of hard thermal loops can 
be absorbed in the variational effective mass. Also, the Goldstone theorem 
is obeyed at any order of the expansion [4]. 

In a recent paper [4], we applied the 2PPI expansion to the calculation 
of the effective potential of A^"^ theory at finite T. At one loop 2PPI, we 
showed equivalence with summation of daisy and superdaisy diagrams [5,6,7]. 
As is well known, at this order of resummation, the phase transition is first 
order. In this paper, we want to do the full two loop calculation of the 
effective potential in the 2PPI expansion. The only two loop 2PPI diagram 
which we have to calculate is the setting sun diagram. As it turns out. 
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this two loop contribution is large enough to turn the first order transition 
into a second order one as expected from lattice simulations and experiment. 
The organisation of this paper is as follows. In section 2, we give a short 
introduction to the 2PPI expansion and its renormalization. In section 3, we 
calculate the one and two loop terms of the 2PPI expansion of the effective 
potential. In section 4, we discuss our numerical results and show that the 
phase transition is second order with mean field critical exponents. 



Let us consider a generic 2PPR diagram (two particle point reducible), i.e. 
a diagram that becomes disconnected when two internal lines meeting at the 
same point are cut (see fig. 1). It is clear that the 2PPR insertions are seagull 
(2 external (f) lines) or bubble graphs. From that, we can conclude that the 
2PPI expansion sums seagull and bubble graphs. How is this summation 
carried out ? We first notice that the seagull and bubble graphs contribute 
to the selfenergy as effective masses and |(0^)c = respectively. 

Therefore, one could naively think that all 2PPR insertions are summed by 
simply deleting the 2PPR graphs from the IPI expansion and introducing 
the effective mass. 



where (p = (0), in the remaining 2PPI graphs. This is not correct since there 
is a double counting problem which can be easily understood in the simple 
case of the 21oop vacuum graph (daisy graph with two petals) of fig. 2. a. 
Each petal can be seen as a self energy insertion in the other, so there is no 
way of distinguishing one or the other as the remaining 2PPI part. We could 
however earmark one of the petals by applying a derivative with respect to 
(p (fig. 2.b). This way the 2PPI remainder (which contains the earmark) is 
uniquely fixed. Now, there are two ways in which the derivative can hit a 
field. It can hit an explicit (p field which is not a wing of a seagull or it can 
hit a wing of a seagull, i.e. an implicit (p field hidden in the effective mass. 
We therefore have 



2 The 2PPI expansion 



m^ = m^ + -(<^^ + A) 



(1) 





^f-(m^^) + A^-^^f-(m^^) 



(2) 
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where T^^^ — S{ip) + Fj^^. Using the definition (1) of the effective mass, we 
can rewrite this as : 

Using an analogous combinatorial argument we have 

and since A/2 = , we find the following gap equation for A : 

(5) 

The gap equation (5) can be used to integrate (3) and we finally obtain : 

T'^\m\ ^) = S{^) + rf ^) - ^ / d^xA' (6) 

with the effective mass given by (1). So indeed, the 2PPR insertions 
can be summed into the effective mass but there is a corection term (the 
last term in (6)) which accounts for double counting. This effective mass 
can be determined self consistently from the gap equation (5) which can be 
rewritten as : 

^ QY2PPI 

^ m'' + -cp'' + X {rff,ip) (7) 

2 om^ 

Up till now, we have remained silent about renormalization. For example, 
A is quadratically divergent and needs renormalization. The question is if af- 
ter renormalization, the simple relation (6) and the gap equation (5) remain 
valid. Fortunately, this is so if one uses a mass independent renormaliza- 
tion scheme. In [2] we showed that using dimensional regularization (which 
eliminates quadratic divergences) and a mass independent renormalization 
scheme, A can be renormalized as : 

An = Z,A + SZ.ip' + (8) 

A 

where renormalization constants come from the standard counterterm La- 
grangian for A0^ : 

«£=f0,*r + f^mV + ^Zx^ (9) 
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and SZ^^^^ contains only the divergences from mass renormalization dia- 
grams which are 2PPR. Using (8), one shows that the gap equation (5) 
becomes after renormahzation : 

with m|, = rf? + |(<^^ + Aj?) and the relation (6) between V"^^^ and V'^^^^ 
becomes 

m') = S{^) + rf/^(mL d^xAl (11) 

The 2PPI effective action is renormalized solely with 2PPI counterterms. 
The other 2PPR counterterms go into the renormalization of the effective 
mass rn^. In minimal subtraction, equation (8) is automatically fulfilled and 
renormalization of the 2PPI expansion becomes extremely simple. 



3 Effective potential at two loop 2PPI 

Prom now on, we will work with the renormalized 2PPI expansion and drop 
the subscript R, keeping in mind that modified minimal subtraction {MS 
scheme) is implicity understood. The first few diagrams in het 2PPI expan- 
sion are displayed in figure 3. We use the imaginary time formalism and 
the Saclay method [8] to calculate the Matsubara sums. For the one loop 
contribution to the 2PPI effective potential, we find : 

V^f^^(m^^) = li:/ln(A;2 + m^) 
whereE/=rE / ^,71' = 47^6" V and 




The effective potential at one loop 2PPI (sum of daisy and superdaisy dia- 
grams) is then given by : 

V;//(m^ if, T) = —if' + -^^ + -—^ [\n ^ - 2 j + QTijff) - -A' (14) 
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where — + ^ (c/?^ + A) is a solution of the one loop gap equation 



—2 2.^f^.^ 

m — m -\ — 

2 2 



(15) 



with 



dm"' J 



i2.r E, 

At two loops, the only 2PPI diagram we have to consider is the setting 
sun diagram. It has been calculated together with the other 2PPR two loop 
diagrams in [9] . For details, we refer to the appendix, especially with respect 
to 2PPI renormalization. For the effective potential at two loop 2PPI, we 
obtain the following expression : 



m 2 -^4 



2^ +4!^ +647r^ 
^A^ - ^^^G.irff) + Gl(m^T) + G,{m\T)){ll) 



where 



3 



2 



2 (47r)^ 



ln^^ + 41n^ + 4+:- 



TT^ 8.966523919\ 



(18) 



+ 



3 p /-^ dq2 

4(27r)^ 



with 



q2nB{q2) , 
12 — = ln| — I 



E. 



?2 



with 



(20) 
(21) 



Y± = {Eg^ + Eg^ + Eq^±g^Y ^ ^ + "^92 "^91 ±92 ) ^ 



(22) 
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The effective mass is a solution of the 2 loop gap equation : 




) 



(23) 



4 Numerical results 

The results obtained in the previous section were used to calculate the ef- 
fective potential for A0^-theory at the 2-loop 2PPI level numerically. These 
calculations involved two steps. The first one is the numerical integration of 
the double integrals in ( |l^) and ( ^). Both integrals should be integrated 
from tot oo. However we found that is was sufficient to approximate the 
infinite integration interval by the compact interval [0,10^]. The stability 
of the numerical results towards changes in this interval and changes in the 
parameters of the model was checked and found to be satisfactory. A second 
difficulty arises from the need to solve the gap equation which is a tran- 
scendental equation. Trying to calculate the effective 2PPI potential as a 
function of 99, we need to solve the gap equation ( ^3]) which gives us the 
variational mass m as a function of 99. This could be done by iteration, but 
this approach will demand a lot of computational ressources. We opted for 
a different approach, which consists of determining as a function of the 
mass, calculate V{(p{m'^),rrf) and plot V{(p) by varying m. 

The parameters in our Lagrangian and the renormalization scale Jt were 
chosen in such a way that comparison with the work of Chiku [13], which 
also uses an effective mass to improve upon the IPI loop expansion, would be 
easy. We refer to this paper for more details about this choice of parameters. 
We used the two parameter sets used by Chiku in [13] and found that the 
results did not depend qualitatively on the particular set. The numerical 
results in this section have been obtained with A = 10.0, = —170 and 



In figure 4 the 1-loop 2PP1 effective potential V{(p) is shown for temper- 
atures in the neighbourhood of the critical temperature. This figure clearly 
indicates a first order transition. Attention should be payed to the fact that 
this first order transition may easily be mistaken for a second order one if 
looked at at a larger scale. Therefore, it is necessary to use a sufficiently fine 



87.6. 
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grid of points around 9^0 where the potential reaches its minimum, and to 
look at temperatures which are sufficiently close to the critical temperature. 
The transition also manifests itself as a discontinuity in the evolution of cpo 
at the critical temperature where ipo jumps from about 0.4 to 0. (figure 5). 

Including the setting sun diagram, the 2PPI effective potential is altered 
to exhibit a second order phase transition, as can be seen from figure 6. 
The second order nature of this transition was extensively checked in the 
neighbourhood of the critical temperature Tc = 24.421. It is also evident from 
the continuous descent of ipo towards zero, where the symmetry is restored 
(figure 7). In figure 7 there seems to be a shoulder structure for which we 
have no physical explanation. 

All of the above mentioned results show qualitative agreement with the 
2-loop analysis of theory in optimised perturbation theory, as done by 
Chiku in [13] 

For the critical exponent /3, which is defined by 

MT) ~ 

for temperatures close to the critical temperature, we found the value /3 = 
0.52 when taking the temperature interval [24, 24.421], indicating consistency 
with the Landau mean-field theory which predicts /3 = 0.5. As a check for this 
last claim, we also investigated the critical behaviour of the second derivative 
of the effective potential and found again accordance with the Landau mean- 
field prediction. 

Summarizing, in this paper we have analysed the \(f)^ theory at finite 
temperature using the 2PPI expansion at the 2 loop order. We found that 
the inclusion of the setting sun diagram which is the only two loop 2PPI 
diagram, changes the phase transition from first order at one loop 2PPI to 
second order at two loop 2PPI. The critical exponents we found were mean 
field. We do not believe that adding more terms of the 2PPI expansion 
will improve the critical exponents. We expect that renormalization group 
resummation of the 2PPI expansion [14] is necessary to breng the critical 
exponents in closer agreement with experiment. 



I' - I' 



Tr 
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Appendix 



In this appendix, we will calculate and renormalize the setting sun diagram. 
Because this is a 2PPI diagram, only the 2PPI counterterms have to be 
taken into account. The 2PPR counterterms are used to renormalize 2PPR 
insertions in 2PPI diagrams. This is a different approach than ordinary IPI 
renormalization where all the 21oop diagrams (two bubble and setting sun 
graphs) have to be considered together and all one loop counterterms have 
to be inserted [10]. In [2] we showed that both approaches are equivalent and 
the purpose of this appendix is mainly to check this explicitely at two loops. 
For those aspects of the calculations which are not related to renormalization, 
we will borrow heavily from [9] to which we refer for further details. 

Using the Saclay method, we find the following expression for the setting 
sun diagram : 

y[Go(m2)+Gl(m^^) + G2(m^^)] (1) 
where Go is the temperature independent part 

r (-^) = f "^""^^ / ^""^^ ^ ^ 1 

"^"^ ^ J (27r)^ J (27r)^ qf + rrf ql + rff {qi + gs)' + rff 

= 3 j d[q,,q2]S{Eg,,Eg,,Ej (2) 

db^^ - (2.)--i(2.)--8E,E,E,3 ^ ^ 

and ^3 = —{qi + and where Gi and G2 are temperature dependent with 
one and two Bose-Einstein factors respectively : 

Gi{m^,T) = 3 J d[qi,q2]nB{qi) [S{E,^, E,,EJ + Si-E,„ E,„ Ej] (5) 



with 



G2{m^,T) = 3 I d[qi,q2]nB{qi)nB{q2)[S{Eq^,Eq2,Eq.^) + S{-Eq^,Eq^,Eq^) 

+ S{Eg^, -Eg^,Eg^) " S {E , E , - E g^)] (G) 
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The setting sun diagram at T = has been calculated in [11,12] in D = 4— 2e : 



Go{m ) = /i 



4e 



{An] 



4-2e ■ 



r(2e - 1) - + 3 - 8.966523919e + ... (7) 



Doing the integral over the angle between qi and q2, G*i(m^,T) can be sepa- 
rated in divergent and finite parts [9] : 



Gl(m^^) 



+ 



(4^ ^ 
3 

4(2^)4 Jo 



2 + In 



dqi- 



E. 
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dq2 

Eq2 



g2ln 



1^ 



-4gi 



with 

and 



X± = {Eq^ + Eq^ + Eq^±q^Y X {-Eq^ + Eq^ + -Egiiga 



.4e 



(27r)3-2^ Eg 



(8) 



(9) 



(10) 



Because 6*2 has a Bose-Einstein factor for each momentumvariable, it is U.V. 
finite and after doing the angular integration one finds : 



G2{m\T) 



4(27r)4 Jq 



dqi- 



E, 



91 







dq2 In 



E 



92 



:ii) 



with 



y± = {Eq^ + Eq^ + Eq^±q^)^ X (--Eg^ + -Ega + Eq^±q^f 
X (Eg^ - Eq^ + Eq^^q^f X (E,^ + Eq^ - -Eg^i^J^ 



(12) 



We now investigate the renormalization of the setting sun diagram at 
finite temperature in some more detail. Since the counterterm Lagrangian is 



(13) 



we find from 



X3\ + crossing 



finite 
finite 



(14) 
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that at one loop 



5Z = 0, 6Zo 



In [2] it was shown that 
so that 



327r 



5Zf ^« = 5Z2 



327r2 



2PPI 



5Z4 — 5Z2 



167r2 



(15) 

(16) 
(17) 



To renormahze the setting sun diagram which is 2PPI, we only have to include 
2PPI counterterms 

O + ^ + ^ = finite (18) 

— 2PPI 2PPI 

This equation determines the two loop 2PPI part of mass renormalization. 
When using the 2PPI counterterm for mass renormalisation of 2PPI dia- 
grams, we should keep in mind to replace the ordinary mass with the 
effective mass m^. This is correct to all orders as was shown in [2]. Coupling 
constant renormalization gives : 



A 

1 A2 



with 



h(rrf) 



2 {An) 



2 

m 

(4^ 

2 

m 



q2 _j_ ^2 



(19) 



(47r)2 

From (A.l) and (A. 7), we have 



In- 



l + 0(e) 



A^ ■m? 



4 (47r) 



+ 



finite 



(20) 



(21) 



As a check on 2PPI renormalization, we find that the ^ In m?/ jj? term coming 
from overlapping divergences gets nicely cancelled and that at two loops 



szl^'P' 



4 {Any 



(22) 
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Analogously, we find from (A.l), (A. 8) and (A. 19) that the temperature 
dependent divergent part gets exactly cancelled by the 2PPI counterterm 
diagram. Doing the renormalization to 0(e*^), one finds that the finite parts 
of Go, Gi and G2 are given respectively by (18), (19) and (21) of section 3. 
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Figure 1: Generic 2PPR diagram 
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Figure 4: The 1-loop effective potential at T=21.5,21.6,21.65,21.7,21.8,22. At 1-loop we 
see a clear first order transition. 



(j)o(T)/(t)o(0) 

1 - 

0.8 
0.6 

0.4- 
0.2- 



Figure 5: The 1-loop calculation of where the potential reaches it's minimum at 
different temperatures shows a discontinuous jump to zero at the critical temperature. 
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Figure 6: The 2-loop effective potential at T = 23.8,24,24.2,24.4,24.6 . Including the 
SS-diagram changes the 1-loop first order transition to a second order one. 
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Figure 7: The 2-loop calculation of 0o shows a continuous approach to zero at the critical 
temperature. 
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